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Stored antiproton beams at the international FAIR facility will provide unique opportunities to 
study hyperons as well as antihyperons in nuclear systems. Precise /-spectroscopy of multi- 
strange hypernuclei will serve as a laboratory for the hyperon-hyperon interaction. Exclusive 
hadron-antihadron pair production close to threshold can measure the potential of a antihadron 
Q relative to that of the coincident hadrons. 

In the present work we explore the production of excited states in double hypernuclei following 
the micro-canonical break-up of an initially excited double hypernucleus which is created by 
the absorption and conversion of a stopped Z~ hyperon. Generally the formation of excited 
hypernuclear states relative to ground states dominates in this model. For different initial target 
nuclei which absorb the S~, different double hypernuclei nuclei dominate. We also compare the 
model predictions with the correlated pion spectra measured by the E906 collaboration. 
In antiproton nucleus reactions the event-by-event transverse momentum correlations of hadron- 
antihadron pairs produced close to threshold contain information on the difference between the 
nuclear potential of the hadron and the associated antihadron. For produced D-meson pairs at 
6.7 GeV/c the sensitivity of the transverse momenta correlation will probably be to small to de- 
duce differences between the potentials for D + and D~ mesons. However, for EE pairs produced 
at 2.9 GeV/c the asymmetry is sufficiently sensitive to predicted differences between the E and 
S potentials even if the momentum and density dependence of the the potential are taken into 
account. 
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1. Introduction 

Quantum Chromo Dynamics (QCD) is the theory of the force responsible for the binding of 
nucleons and nuclei and thus of a significant fraction of the ordinary matter in our universe. While 
the internal structure of hadrons and the spectra of their excited states are important aspects of 
QCD, it is at least equally important to understand how nuclear physics emerges in a more rigorous 
way out of QCD and how nuclear structures - nuclei on the small scale and dense stellar objects 
on the large scale - are formed. In particular the role of strangeness in neutron stars has not been 
settled yet. 

The study of hypernuclei can illuminate features that are obscured in conventional nuclear 
systems. The hyperon offers a selective probe of the hadronic many-body problem as it is not 
restricted by the Pauli principle. Since direct scattering experiments between two hyperons are 
impractical, the precise spectroscopy of multi-strange hypernuclei provides a unique chance to 
explore the hyperon-hyperon interaction. Significant progress in nuclear structure calculations in 
chiral effective field theory nurtures the hope that detailed information on excitation spectra of 
double hypernuclei will provide unique information on the hyperon-hyperon interactions. At the 
same time, the A-N and A-A weak interaction can be studied by hypernuclei decays, opening a 
inimitable window for the four-baryon, strangeness non-conserving weak interaction. 

Furthermore, the physics of strangeness in hadronic systems is a constantly developing field. It 
brings up new, often unexpected results, new challenges and open questions like anti-hypernuclei, 
kaonic nuclei, exotic hadronic states like the controversially discussed pentaquark baryons and the 
H-dibaryon. Based on G-parity transformation [1] Diirr and Teller predicted within an early form 
of a relativistic field theory a strongly attractive potential for antiprotons in nuclei [2]. First in- 
vestigations of antiproton-nucleus scattering cross sections [3, 4, 5] showed however disagreement 
with a strong attractive potential. Later, X-ray transitions in antiprotonic atoms [6, 7, 8, 9, 10] 
gave also hints for an attractive potential albeit with large uncertainties [11, 12]. More comprehen- 
sive studies [13] of antiprotonic X-rays as well as recent analyses of the production of antiprotons 
in reactions with heavy ions resulted in attractive real potentials in the range of about -100 to - 
150 MeV [14, 15, 16]. It is obvious that G-parity can only provide a link between the AW and NN 
interactions for distances where meson exchange is a valid concept [17, 18]. For distances lower 
than about 1 fm, quark degrees of freedom may play a decisive role. Nevertheless this considera- 
tion still indicates that a direct comparison of the interactions of baryons and antibaryons in nuclei 
may help to shed some light on the nature of short-range baryon-baryon forces. 

2. Production of Double Hypernuclei 

The simultaneous production and implementation of two A particles into a nucleus is intricate. 
There is a possibility to produce multi-strange hypernuclei in heavy ion collisions via coalescence 
[19, 20]. The first observation of antihypernuclei by the STAR collaboration impressively illustrates 
the potential of this method [21]. However, high resolution spectroscopy of excited states is not 
feasible. To produce double hypernuclei in a more 'controlled' way the conversion of a captured 
S~ and a proton into two A particles can be used. This process releases - ignoring binding energy 
effects - only 28 MeV. For light nuclei there exists therefore a significant probability of the order 
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Figure 1: Present knowledge on S=-l nuclei (blue) and S=-2 nuclei (yellow). Only very few individual 
events of double hypernuclei have been detected and identified so far. First antihypertritons were recently 
observed by the STAR collaboration [21]. 

of a few percent that both A hyperons are trapped in the same excited nuclear fragment [22, 23, 24, 
25, 26, 27]. 

Unfortunately S~ hyperons produced in reactions with stable hadron beams have usually 
rather high momenta. Therefore, direct capture of the S~ in the nucleus is rather unlikely. Even 
in case of the (K~,K + ) double strangeness exchange reaction, S~ hyperons are produced with 
typical momenta of 500 MeV/c at beam momenta around 1.8 GeV/c [24, 28]. The advantage of 
this production process is that the outgoing K + can be used as a tag for the reaction. A drawback 
is the low kaon beam intensity and hence the need for thick primary targets. Furthermore, as a 
consequence of the large momentum transfer, the probability to form bound H~ states directly is 
rather small on the level of 1% [29, 30] and the production of quasi-free S~ dominates. Still the S~ 
hyperons in the quasi-free region may be absorbed into the target nucleus via a rescattering process 
on a nucleon which itself is knocked out of the primary nucleus. This two-step process is predicted 
to exceed the direct capture by more than a factor of 6 [23]. 

On the other hand most (~ 80%) S~ hyperons escape from the primary target nucleus in 
(K~ , K + ) reactions. However, in a second step, these S~ hyperons can be slowed down in a dense, 
solid material (e.g. a nuclear emulsion) and form S~ atoms [31]. After an atomic cascade, the 
S-hyperon is eventually captured by a secondary target nucleus and converted via the Z~p —> AA 
reaction into two A hyperons. In a similar two-step process relatively low momentum S~ can 
also be produced using antiproton beams in pp —> S~S or pn —> S~H° reactions if this reactions 
happens in a complex nucleus where the produced S~ can re-scatter [32, 33]. The advantage as 
compared to the kaon induced reaction is that antiprotons are stable and can be retained in a storage 
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ring. This allows a rather high luminosity even with very thin primary targets. 

Because of the two-step mechanism, spectroscopic studies, based on two-body kinematics 
like in single hypernucleus production, cannot be performed. Spectroscopic information on double 
hypernuclei can therefore only be obtained via their decay products. The kinetic energies of weak 
decay products are sensitive to the binding energies of the two A hyperons. While the double 
pionic decay of light double hypernuclei can be used as an effective experimental filter to reduce 
the background [34] the unique identification of hypernuclei groundstates only via their pionic 
decay is usually hampered by the limited resolution. Instead, y-rays emitted via the sequential 
decay of excited double hypernuclei may provide precise information on the level structure. 

2.1 Statistical Decay of excited Doubly Strange Nuclei 

The PANDA experiment [33] which is planned at the international Facility for Antiproton and 
Ion Research FAIR in Darmstadt aims at the high resolution y-ray spectroscopy of double hyper- 
nuclei [32]. An important question is to what extent double hypernuclei in excited, particle stable 
states are populated following the break-up of an highly excited doubly strange nucleus which is 
formed after the absorption and conversion of a stopped S~. For light nuclei even a relatively 
small excitation energy may be comparable with their binding energy. We therefore consider the 
explosive decay of the excited nucleus into several smaller clusters as the principal mechanism of 
de-excitation. 

To describe this break-up process we have developed a model [35] which is similar to the 
famous Fermi model for particle production in nuclear reactions [36]. We assume that the nu- 
cleus with mass numbers Ao, charge Zo, and the number of A hyperons Ho (here Ho=2) breaks-up 
simultaneously into cold and slightly excited fragments, which have a lifetime longer than the de- 
cay time, estimated as an order of 100-300 fm/c. In the model we consider all possible break-up 
channels, which satisfy the mass number, hyperon number (i.e. strangeness), charge, energy and 
momentum conservations, and take into account the competition between these channels. 

The excitation energy of the initial, highly excited double A nucleus is determined by the bind- 
ing energy of the captured S~ hyperon. Unfortunately, still very little is established experimentally 
on the interaction of H~ hyperons with nuclei. Various data suggest a nuclear potential well depth 
around 20MeV (see e.g. [37, 10]). Calculations of light 2 atoms [31] predict that the conversion 
of the captured S~ from excited states with correspondingly small binding energies dominates. In 
a nuclear emulsion experiment aS" capture at rest with two single hyperfragments has been ob- 
served [38] which was interpreted as S~ + 12 C + ^Be reaction. The deduced binding energy 
of the S~ varied between 0.62 MeV and 3.70 Me V, depending whether only one out of the two 
hyperfragments or both fragments were produced in an excited particle stable state. In order to 
take into account the uncertainties of the excitation energy of the converted S~ -states, the calcula- 
tions were performed for a range of energies < Bz < E max , constructing in this way the excitation 
functions for the production of hypernuclei. 

Fig. 2 shows as an example the production of ground (g.s.) and excited (ex.s.) states of single 
+ one free A (SHP), twin (THP) and double (DHP) hypernuclei in case of a 12 C target as a function 
of the assumed S~ binding energy. With increasing 3~ binding energy the excitation energy of the 
excited primary ^\B* nucleus decreases from left to right from about 40 MeV to 15 MeV. For all 
excitation energies above 20 MeV the production of excited double hypernuclei dominates (green 
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Figure 2: Predicted production probability of ground (g.s.) and excited states (ex.s.) in one single (SHP), 
twin (THP) and double hypernuclei (DHP) after the capture of a Z~ in a 12 C nucleus and its conversion into 
two A hyperons . The lower and upper scale shows the binding energy of the captured Z~ and the excitation 
energy of the initial J^B nucleus, respectively (from Ref. [35]). 



triangles). This can be traced back to the opening of several thresholds for various excited double 
hypernuclei already at moderate excitation energies. Only for small binding energies and hence 
large excitation energies the production of single and twin hypernuclei is significant (~10%). The 
non-monotonic behaviour for single hypernucleus + one free A production reflects the fact that 
the various lowest thresholds are relative high and widely separated, e.g. ^ 2 B+A at Bs =23. 9 MeV 
followed by J^B+n+A at 11.3 MeV. Twin-hypernuclei are only produced for S~ binding energies 
below the threshold for ^Li + ^He with Bs=13.6 MeV. As discussed above the frequent observation 
of twin-hypernuclei [39, 40, 41, 42, 38, 43] signals a conversion from a S state with only moderate 
binding energy. In this range of Bs the production probability of double hypernuclei is comparable 
to previous estimates within a canonical statistical model [23, 24]. It is important to stress that these 
numbers do not include a possible pre-equilibrium emission of hyperons during the capture and 
conversion stage. With respect to the number of stopped S~ hyperons, pre-equilibrium processes 
will decrease the yield of double hypernuclei relative to the yield for single hypernuclei (see e.g. 
[27]). Indeed in the simulations for the planned PANDA experiment [34] a joint capture x conversion 
probability of 5% was assumed to mimic this pre-equilibrium stage. 

Using different S-absorbing stable secondary target nuclei 9 Be, 10 B, n B, 12 C and 13 C one finds 
that different double hypernuclei dominate for each target [35]. Thus combining the information 
shown in Fig. 2 with the measurement of two pion momenta from the subsequent weak decays a 
unique assignment of various newly observed y-transitions to specific double hypernuclei seems 
possible as intended by the PANDA collaboration [32, 34]. 
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2.2 The E906 Puzzle 

In 2001 the BNL experiment E906 reported the observation of the AA H hypernucleus by mea- 
suring the sequential pionic decays after a (K~,K + ) reaction deposited two units of strangeness in a 
9 Be target [44] (see Fig. 3). Two structures in the correlated 7i~ momenta at (133,114) MeV/c and 
at (1 14,104) MeV/c were observed. The first structure was interpreted as the production of A H+ A H 
twins while the bump at (114,104) MeV/c was attributed to pionic decays of the double hypernu- 
cleus AA H. However, as it was pointed out by Kumagai-Fuse and Okabe also twin A-hypernuclear 
decays of A H and A He are a possible candidate to form this peak if excited resonance states of 6 Li 
are considered [45]. More recently Randeniya and Hungerford showed that the published E906 
data can be reproduced without the inclusion of decay and that it is more likely that the decay 
of J^He was observed in the E906 experiment [46]. In their analysis this double hypernucleus was 
accompanied by a background of coincident decays of single hypernuclei pairs A H+ A H, A H+ A H, 
and A H+ A H, respectively. 

Fig. 4 shows the predicted relative probabilities for the production of particle stable twin 
and double hypernuclei in the E906 experiment after the capture and conversion of a stopped S~ 
in a secondary 9 Be target S~+ 9 Be— ^Li*. As before a S~ binding energy of 0.5 MeV seems 
reasonable corresponding to an AA Li excitation of about 29 MeV. The produced yields for this 
excitation energy are shown in Fig. 3. Here, ground state and excited state(s) - if they exist - have 
been added and the pion momenta for groundstate decays are assumed. Note, that the production 
of A H+ A H twins is even at an excitation energy of 35 MeV energetically not possible and - unlike 
to a canonical calculations [24] - does therefore not occur in our micro-canonical model. 

Let us first discuss the structure at (114,104) MeV/c which has been attributed to double hy- 
pernuclei decays. Generally the production of double hypernuclei is energetically favored over the 
production of twins: all possible channels with twin-production lie energetically significantly above 
the thresholds for AA Li, J^He and ^He production in case of the ]^Li compound picture. Cor- 
respondingly the production of A H and A He twins which has been suggested as a possible source 
of the peak structure around (114,104) MeV/c [45], is in our model by a factor of 30 lower than 
the j^He production probability. Unlike in the canonical model of Ref . [24] the J^He probability 
exceeds also the one of a ^H by more than two orders of magnitude and the production of ^H 
by a factor of about 17 in our model. Of course, a direct comparison with the E906 data requires 
a detailed consideration of the branching ratios for pionic two-body decays many of which are 
not known so far. Keeping that caveat in mind our microcanonical model supports independently 
on the assumed production scheme the interpretation of the E906 observation by Randeniya and 
Hungerford [46] in terms of J^He decays. Decays from the ground or even exited states of AA Li or 
AA Li can possibly contribute some background to the structure at (1 14,104) MeV/c. 

In order to describe the (133,114) MeV/c structure of the E906 experiment, the production 
of A H+ A H twins seems mandatory [44]. However, the A H+ A H+t mass lies above the initial mass 
mo = m(H~) +m( 9 Be) and can therefore not be produced in the S~+ 9 Be compound production 
scheme (see right part of Fig. 4). With an energy of 12.6 MeV below mo, the channel A H+ A He is 
the most likely twin in the present scenario, followed by the A H+ A He+n decay. Given however 
the experimental precision of about 1 MeV/c for the momentum calibration in E906 [44], neither 
the decays of A H+^He with (133,108) MeV/c nor decays of A H+^ He pairs with (133,99) MeV/c 
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Figure 3: Momenta of two correlated pions measured by the E906 collaboration (grey squares). Overlaid 
are production yields of various twin pairs (blue), double hypernuclei including their excited states (green) 
and double hypernuclei with a stable groundstate only (red). The areas of the circles are proportional to the 
production yields predicted by the statistical model. Pionic decay probabilities are not included in this plot. 



seem to explain the structure around (133,114) MeV/c. Of course in particular the first one could 
contribute to this enhancement in the tail region. Even more intriguing is however the fact that 
the statistical models predict a production of A H+ A H twins exceeding the A H+ A H production. 
Considering furthermore the branching ratios for two-body 7i~ decays of ^ 7t -^ He /T tota i ^0.26 
[47] and T n - + A He /T tota i « 0.5 [48] the absence of a bump which could be attributed to A H+ A H 
is particularly puzzling. A similar conclusion is obtained [35] within the quasi-free/rescattering 
picture of Yamamoto et al. [24] resulting in the production of excited ^He or AA H nuclei. 

At first sight it seems that an alternative production process than the ones discussed so far 
is required to explain the singular structure at (133,114) MeV/c in terms of A H+ A H twins. Note 
however that in the initial analysis of the E906 data the bump at (133,114) MeV/c served as a 
calibration point for the pion momenta, taking the decay of A H+ A H twins as granted [44, 49]. If 
that structure were indeed caused by A H+ A He twins with (133,108) MeV/c, it would of course 
influence the momentum scale in the region of the (1 14,104) MeV/c bump. This bump would then 
be shifted to approximately (108,97) MeV/c. Considering the uncertainty of AZ?aa also such a 
momentum scale would be compatible with the decay of AA He and ^Li or a mixture of both. The 
absence of AA Li nuclei may be related to the decreasing pionic decay probability with increasing 
nuclear mass. Clearly, the present statistical decay model needs to be complemented by quantitative 
weak decay calculations (see e.g. Ref. [50]) to further corroborate our conjecture. 
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Figure 4: Relative production probability of double (top part) and twin hypernuclei (lower part) for a pri- 
mary J^Li nucleus as a function of its excitation energy. 



3. Antihyperons in Nuclei 

Concerning antibaryons, realiable information on their nuclear potential are available only for 
antiprotons. Antihyperons annihilate quickly in normal nuclei and spectroscopic information is 
therefore not directly accessible. Mishustin recently suggested to study deeply bound antibaryonic 
nuclei via various characteristic signals in their decay process [51, 52]. It is however not obvious 
whether these proposed observables will provide unique and quantitative signals of deeply bound 
antihyperonic systems. 

Quantitative information on the antihyperon potentials relative to that of the corresponding 
hyperon may be obtained via exclusive antihyperon-hyperon pair production close to threshold in 
antiproton- nucleus interactions [53]. Once these hyperons leave the nucleus and are detected, their 
asymptotic momentum distributions will reflect the depth of the respective potentials. However, 
since in the pp center-of-mass the distribution of the produced baryon-antibaryon pair will usually 
not be isotropic, the analysis can rely only on the transverse momenta of the outgoing baryons. In 
Ref. [53] it was demonstrated that the individual transverse momentum distributions alone do not 
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allow to extract unambiguous information on the potential of antihyperons. Instead, the transverse 
momentum asymmetry aj defined in terms of the transverse momenta of the coincident particles 
can be used to explore event-by-event correlations: 

aT = PT(A)- PT (K) ^ (31) 
p T (A)+p T (A) 

In a purely classical, non-relativistic picture this asymmetry is of the order of AU /A- Eo, where AU 
is the potential difference and Eq the typical kinetic energy of the hadrons. 

3.1 Antihyperon-Hyperon Pair Production and Propagation in Nuclei 

In Ref. [53] we have examined the influence of the potentials on the transverse momentum 
asymmetry of coincident hyperons and antihyperons by means of a schematic Monte Carlo simu- 
lation. Albeit crude, this classical approach allows to explore the role of different features of the 
reaction in a transparent way. 

The absorption of the antiprotons entering the target nucleus determines the points of annihi- 
lation inside the nucleus and the paths which the eventually produced hyperons and antihyperons 
have to pass inside the nucleus prior to emission. Because of the strong absorption of the antihyper- 
ons, the em//?edhyperon-antihyperon pairs are - unlike in inclusive reactions [54, 55] - created close 
to the corona of the target nucleus at a typical density of 20 to 25% of the central nuclear density. 
In reactions close to threshold the Fermi motion of the protons inside the nuclear target contributes 
significantly to the final momenta. Lacking any detailed experimental information it is assumed 
that the annihilation cross sections for antihyperons show a similar momentum dependence as the 
pp system [5, 56]. 

The energy and the momentum of the baryons propagating within the nucleus are related 
according to [57]: 

(£-V) 2 = (M + S) 2 + P in 2 (3.2) 

Here V and S denote the real part of the vector and scalar potential, respectively. The relation 
between the momenta inside and outside of the nuclear potential are approximated by 



p 2 

"out 



+ M 2 = (^(M + S) 2 + P in 2 + V) 2 . (3.3) 



Refractive effects at the potential boundary were ignored. The default parameters for the scalar and 
vector potentials of the various baryons at normal nuclear density po were adopted from Ref. [58]. 
Since the antiproton annihilation and the subsequent antihadron-hadron pair production take place 
in the nuclear periphery at low densities p we assumed a linear density dependence °c p /p for all 
vector and scalar potentials. 

3.2 Transverse Momentum Correlations at PANDA 

For a compact representation of event-by-event correlations we examine the transverse mo- 
mentum asymmetry OLj as a function of the longitudinal asymmetry Gfc, where Gfc is defined by 
analogy to Eq. 3.1 in terms of the longitudinal momentum components: 

^= PL \T PLi Sv 

p L (A)+p L (A) 
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Figure 5: Left part: Average transverse momentum asymmetry of AA(solid line), SS (dashed line) and 
D + D~ pairs (dotted line) produced exclusively in 1.66, 2.9 and 6.7 GeV/c p 12 C interactions, respectively. 
Right part: the same as in part (a) but taking the momentum dependence of the potential into account. 

In future more sophisticated analysis methods which e.g. use mixed events as a reference, further 
details of these correlations may be revealed. 

Transverse momentum correlations can in principle be analyzed for each hadron-antihadron 
pair produced exclusively in pA interactions. As examples, the solid, dashed and dotted histograms 
in figure 5a show the average transverse momentum asymmetries of AA, SS and D + D~ pairs 
produced in 1.66, 2.9 and 6.7 GeV/c p + 12 C interactions, respectively. For simplicity, isotropic 
center-of-mass distributions were assumed in case of the SS and D + D~ production. For the S and 
S baryons the same absorption cross sections as for the A and A were adopted, whereas for D~ 
and D + mesons energy independent absorption cross sections of 10 and 90 mb, respectively, were 
taken. In the left part of Fig. 5 no momentum dependence of the potentials was considered. 

In line with the classical picture mentioned above {aj AU /4 • Eq) the smaller asymmetries 
for the heavier particles are a consequence of the the large S hyperon and D meson laboratory 
momenta. In the simulations of Fig. 5b the momentum dependence of the potentials was added by 
means of the scaling factor [59]. As expected the transverse asymmetry is shifted towards more 
negative values. Nevertheless the asymmetry remains sizable even for the SS case. At PANDA 
when reaching the design luminosity a measurement of a? for S3 pairs with a precision of about 
10% will require a measurement time of typically less than a day. Furthermore it puts only moderate 
constraints on the detector performance, e.g. the tracking capabilities of the central vertex detector. 

The fact that energy and momentum conservation are the main ingredients of the proposed 
method raises hope that similar results might be obtained by more sophisticated calculations. Since 
most of the emitted hyperon-antihyperon pairs are created in the nuclear periphery at subsaturation 
density, a neutron skin of neutron rich target nuclei may help to explore different effective poten- 
tials. Possible deflections at the potential boundary which are ignored in the present work may be 
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at least partly eliminated by demanding that the target nucleus remains intact. Furthermore, it may 
be interesting to study questions related to the formation time by using target nuclei of different 
size. Finally, it should be noted that this method can be extended to the case of photo- or even 
electro-production of short-lived resonances in nuclei decaying into particle-antiparticle pairs. 
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